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ABSTRACT 

The behavior of the system hazard rate function of a two 
component parallel system is investigated. The inter- 
relationships between the probabilities that the components 
composing the system are alive and the system hazard rate is 
examined with special attention to certain points where 
there are important changes in the behavior of the hazard 
rate function. The behavior of the system hazard rate 
function is shown to depend upon the rates of change of the 
probabilities that the components of the system are alive. 
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I. INTRODUCTION 



In the field of reliability, there has been an extensive 
study of the properties of the increasing hazard rate (IHR) 
and decreasing hazard rate (DHR) classes of life distribu- 
tions. The exponential life distributions with their 
constant hazard rates form a dividing line between these two 
classes. It is reasonable to expect that the life distribu- 
tion of redundant systems composed of devices with exponen- 
tial life distributions fall into the IHR class. However, 
when two devices with different exponential life distributions 
are combined in a parallel redundant system, the system 
hazard rate is at first increasing and then decreasing. The 
underlying causes of this interesting behavior is the subject 
of the thesis. 

The usual definition of the probabilistic hazard rate 
function of a device with a random lifetime T will be used 
throughout this thesis. The hazard rate r(t) of a device 
is the density function f(t) for the lifetime of the device 
divided by its survival function F(t)=P[T>t]. 

A distribution has an increasing hazard rate if r(t) is 
monotone non-decreasing in t, and has a decreasing hazard 
rate if r(t) is monotone non-increasing in t. References 1 
and 2 contain an extensive survey of the properties of the 
IHR and DHR classes of life distributions. The exponential 
life distributions are included in both of these classes. 
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In Reference 3, Esary and Proschan gave sufficient conditions 
for a system to have an increasing hazard rate when it is 
composed of identical components with increasing hazard 
rates. A parallel system of two components, with identical 
exponential life distributions, satisfies these conditions 
and it has an increasing hazard rate. Esary and Proschan, 
also gave an example to show that when two components with 
different exponential life distributions form a parallel 
system, the system survival function need not have an 
increasing hazard rate. In this case the system hazard rate 
initially increases, overshoots and then decreases to an 
asymptotic value equal to the lowest of the two component 
failure rates. 

In Section II, the derivation of the system hazard rate 
for a system of two parallel components having constant 
non-identical failure rates is reviewed. In Section III, 
a characterization of the system hazard rate as a function 
of some conditional state probabilities is presented and 
their interrelationships are discussed. In Section IV, an 
empirical approximation for the time at which the system 
hazard rate is a maximum is exhibited. 
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II. SYSTEM HAZARD RATE 



Consider a system of two independently functioning 
components with exponential lifetimes T^,T 2 and survival 
functions 

-X,t 

F 1 (t ) = e 

-X ? t 

F 2 (t) = e 

Assume without loss of generality that X 2 >X^. The equations 
for the system survival function and the system hazard rate 
are generally known. 

With the components in parallel, the system lifetime T 
has a survival function 

F(t) = P[max T 1 ,T2 > t] 

-x 1 t -x 2 t -(x 1 +x 2 )t 

= e +e • -e 



The density of T is 

•fCt) = - P(t) 

-X.t -X„t -(X,+Xp)t 

= X x e +X 2 e -(X 1 +X 2 )e 



The system hazard rate is 



r(t) 



f(t) 

F(t) 

-X,t -X„t -(X 1 +X,)t 

X 1 e +X 2 e d -(X 1 +X 2 )e 1 



F(t ) 
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Examining r(t) at t=0 and as t approaches infinity gives 



the expected results; 

r(0) = 0 

lim r(t) = min(X,,X ? ) = X 
t-*» 1 ^ 



The derivative of r(t) with respect to time is 



r ’ (t ) = 



F(t)f'(t) + f^(t) 



F 2 (t) 



o - X, t n — X 0 1 n — (X-.+X_)t 

[X 2 d e 1 tX^e 2 -(X 1 -X 2 )^]e 1 2 

F 2 (t) 



Looking at the sign of r'(t), it is clear that the 

hazard rate is increasing if and only if 
2 ~X t p — Xp t 2 

X 2 e +X^ e > (X -X ) and decreasing if and only if 

2 “ X-i t 2 — Xp t 2 

X 2 e +X 1 e < ( X i“ X 2 ) * If X 1 =X 2 5 r(t) ls 



increasing for all t>0. The system hazard rate must be at 
a maximum for the value of t which is a solution to 



x 2 Xlt +x 

A 2 e +A^ 



2 e ^-(X ^) 2 



0 
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III. SYSTEM HAZARD RATE 
IN TERMS OP 
STATE PROBABILITIES 



If the parallel system is alive at some time t>0, it 
exists in one of three possible states, i.e. both components 
alive, component one alive and component two dead, or 
component one dead and two alive. The system hazard rate 
and its derivative have simple expressions in terms of the 
conditional probabilities for these states and their 
derivatives. The probability that the system is in state 
i, i=0,l,2, given it is alive at a time t, will be called 
the state probability, denoted by P^(t). For a system with 
lifetime T and component lives T^ and T 2 , the state 
probabilities are: 

(1) both components alive 

P 0 (t) = P[T 1 >t ,T 2 >t |T>t] 

P [ T^ >t , T 2 >t , max ( T^ , >t ) ] 

~~ PCmaxCT^ ,T' 2 >t ) J 

-(x 1 +x 2 )t 

e 

F(t) 

(2) component 1 alive and 2 dead 

P x (t) = P[T 1 >t,T 2 <t |T>t] 

P[T. >t ,T 2 <t jmaxCT^^ ,T 2 >t ) ] 

P[max(T-^ ,T 2 >t ) ] 
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P 1 (t) = 



_X l t -( x 1 +x 2 ^ t 
e -e 

F(t) 



= 1 - 



-x 2 t 



F(t ) 



X_t 

(e 2 -l)P n (t) , 



(3) component 2 alive and 1 dead 



P 2 (t) = P[T 1 <t ,T 2 >t |T>t] 



P[T 1 <t ,T 2 >t jinaxCT^ ,T 2 >t ) ] 
PCmaxCT^ 3 T 2 >t ) ] 

-x 2 t -(x 1 +x 2 )t 

e -e 



= 1 - 



F(t ) 



-x 1 t 



F(t) 



X. t 

(e 1 -l)P Q (t) . 



It is clear that Pg(t) + P^Ct) + Pp(t) = 1, for all t. 

Examining the state probabilities at t = 0 and as t 
approaches infinity produces the following results. 



V 0) 



P 1 (0) 



P 2 (0) 



lj lim Pq ( t ) 

t->-°° 

0 3 lim P, (t ) 
t-*» 1 

0, lim Pp(t) 

t->co d 



= 0 



= 1 



= 0 . 
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At first glance, it might seem that P^Ct) would equal 
PgCt) at the me'an Xxifife' to failure for component 2, hut, 
p i(t) = P Q (f) if and only if 



~x t -(x 1 +x 2 )t 

e -e 



F(t ) 



-(x 1 +x 2 )t 

F(t) 



or 



or 



or 



Similarly, 



P 2 (t) 



-x 1 t -(x 1 +x 2 )t 

e = 2e 



1 = 2e 



-x 2 t 



t = t— In 2 



P Q (t) if and only if 



t * t— In 2 



The probability that a particular single component is alive 
at some time t is greater than the probability that both are 
still alive for any time t greater than the mean tine to 
failure of the other component multiplied by a constant In 2. 

Further, for X 2 >X^, P^(t) > P 2 (t) for all t, since 

P-^(t) > P 2 (t) if and only if 

-x 2 t -x 1 t 

1 - > 1 - 

F(t ) F(t) 
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or 



-Apt -At 
e < e 1 



> 



or 



X 



2 



> X 



1 * 



The derivatives of the state probabilities with respect 
to time are: 



P 0 (t) = 



-(X,+X ? )t -( X-i +Xp )t 

P(t)[-(X 1 +X 2 )e -e P (t) 



F 2 (t) 



-(X, +X ? )t 

[ '^1 - (X 1 +X 2 ) ] 



p(t) 



F(t) 



= Pq ( t ) [r ( t ) - (X 1+ X 2 )] . 



P 1 (t) =- 



-X t -X t 

P(t)(-X 2 e ) - e 2 F (t) 

F 2 (t) 



-x 2 t 



F(t) 



[X 2 - r (t ) ] . 



P 2 (t) = - 



A - 1 ^ A ^ t | 

F(t)(-X ie 1 ) - e 1 F (t) 

F 2 (t ) 






F(t) 



[X, - r (t ) ] . 
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P n (t) decreases in time from P Q (0) = 1 to lim P q ( t ) = 0, and 
, t-+°° 

pQ(t)<0 for all t. P 1 (t ) is increasing in time from 

P,(0) = 0 to lim P 1 (t ) = 1, and P-^(t)>0 for all t. 
t^°° 

P 2 (t) has a maximum at a time x which can be found by 

t 

substituting for r(t) in P 2 (t) and setting this expression 
equal to zero, i.e. when 



-(X-i+Aq 

(X 1 -X 2 )e 



Multiplying both sides by e 
gives 



)t -(2X 1 +X 2 )t 

+ X 2 e 

P 2 (t) 



0 . 



(x 1 +x 2 )t 

and combining terms 



“A-^t X 2 •* X i 



and 



X 1 ^ (x 2" X l ) * 



The value of P 2 (t) at maximum is 



1 + — 

x 2 2 - x 1 x 2 - (x 2 -x 1 ) 1 

x 2 2 - X 1 X 2 + X 2 X 2 //X 1(X 2 -X 1 ) - (A 2 -A 1 ) 1+(X 2 //X 1 ) 



and this is less than one for all ^2 >X l * 

As a consequence of the definition of the state 
probabilities and the requirement that 

P 0 (t) + P x (t) + P 2 (t) = 1 , 
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the derivatives are related and 



Pg(t) + P|(t) + PgCt) = 0 . 

Therefore, P ^ t ) = -[Pg(t) + P^(t)] and P 2 (t) increases 

! I t t 

when -P Q (t ) > P^(t) and it decreases when -Pg(t) < P^(t) . 
P 2 (t) increases as long as Pg(t) decreases in time faster 

than P 1 (t) increases, reaches its maximum value when 

1 1 

-Pg(t) = P^(t), and then decreases because P^(t) is increas- 
ing in time faster than Pg(t) is decreasing. 

A hazard rate may be thought of as a conditional 
instantaneous probability of failure at time t, given 
survival to time t. The hazard rate for the system, given 
that it is in one of its three states, is 0, A^, and X 2 
respectively. The system hazard rate may be written in 
terms of these hazard rates and the state probabilities as 

r (t ) = 0-P 0 (t) + X 1 P 1 (t) + X 2 P 2 (t) 

-(x 1 +x 2 )t — x_t -(x 1 +x 2 )t -Xgt -(\ 1 +x 2 )t 

= 0 • + A. ^ — + X_ 

F(t ) 1 F(t) d F(t) 

which reduces to the expression previously given for r(t). 

The derivative of the system hazard rate with respect to 
time may be written as 

r (t) = X 1 P 1 (t) + X 2 P 2 (t) 



1H 



Then the condition for the maximum point of the hazard rate, 
r(t) = 0, requires 



A-jP^Ct) + ( t ) = 0 • 

f t 

Substituting for P^(t) and P2(t) yields 

-X 2 t -A^ 

A-, [A, - r (t ) ] + X 0 - [ A-. - r(t ) ] = 0 . 

1 F(t) 2 2 P(t) 1 



Solving this for r(t) gives its maximum value 



r(t) = 



-A,t -A 9 t 

X l X 2*- e 1 + e J 

-At 

X l e + X 2 e 1 



Prom P 1 (t ) = 



-A 2 t 



F(t) 



[A 2 - r(t)] > 0 for all t, it is cleai 



that r(t) < A 2 for. all t. At its maximum value r(t) > A- 
if and only if 



-At -Apt 

A-, Ap [e 1 + e 2 ] 

— rr. - > \ 

-At -A,t A 1 

X l e ^ + X 2 e 



or 



-A, t -Apt p -Apt -A, t 

X l X 2 e +A 1 A 2 e > A x e 2 +x ± x i e 1 > 
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or 



A^ ^2 ^ A 

or A 2 > A 

Therefore, at its maximum. 

Prom p'(t) = 3 [A, - 

^ F(t) 1 

r(t) = A^ at the time t = 

its maximum value. If t > t 
and only if • 



2 

1 5 
1 * 

A 1 < r(t) < A 2 

r(t)] , it is clear that 

Y~ In ( r -| r -) when P 2 (t) 

, then r(t) > A-^. r(t) > 



A 1 P 1 (t) + A 2 P 2 (t) > X 1 . 

or A 2 P 2 (t) > A 1 (l - P 1 (t)) , 

or (A 2 -A 1 )P 2 (t) > A 1 P Q (t) . 

Since 

At 

P 2 (t) = (e 1 - 1 ) P Q (t) , 

then r(t) > A^ if and only if 

A, t 

(A 2 -A 1 )(e 1 - 1 ) > A x , 

and for t > t = In ( •. — ) , 

A 1 A 2 A 1 



reaches 
A x if 
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it follows that 




ly- ln(X 2 -X 1 )] 



- 1 ) 




X 



1 



As X 2 approaches X-^ in value x becomes very large, and if 
X-^ = X 2 = X , the system hazard rate r(t) approaches X 
asymptotically from below. For X 2 >> X-^ , x becomes very 
small and the system hazard rate is greater than X^ almost 
instantly . 

From r’(t) = X^P^(t) + X 2 P 2 (t) , r(t) is increasing if 

f f 

and only if X^P^(t) + X 2 P 2 (t) > 0 , and decreasing if and 

t f 

only if X-^P^(t) + X 2 P 2 (t) < 0 . The behavior of the system 
hazard rate is directly dependent upon the rates of change 
of the state probabilities P-^(t) and P 2 (t). r(t) increases 
as long as P^(t) and P 2 (t) increase and it continues to 
increase as long as X^P^Ct) > I ^ 2 P 2 ( t ) I • F° r 
XiPi(t) < | X 2 ? 2 ( t ) | , r(t) decreases. 

Figure 1 shows the state probabilities and the system 
hazard rate plotted as a function of time. As the failure 
rate of component two increases in relation to the failure 
rate of component one, these curves are ’’pushed towards the 
left", the system hazard rate increases rapidly, and its 
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maximum value increases. PgCt) decreases to zero very 
rapidly while P-^(t) increases to one almost instantly. 
P£(t) becomes more peaked and it reaches its maximum value 



sooner. 



r(t) 




Lf\ 



II 

CVJ 



II 

T— 

r< 



00 
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FIGURE 1. r(t), P, (t), and P ? (t) versus time. 



IV. EMPIRICAL RELATIONSHIP 



BETWEEN TWO MAXIMA 



An empirical investigation of the relationships between 
the state probabilities and the system hazard rate was also 
undertaken. The system equations were solved for various 
values of their parameters which permitted comparison of 
the resulting curves and verification of the analytic 
results of Section III. 

The time at which the system hazard rate reaches its 

maximum, T , can be found by solving 
s 



x 



2 

2 




2 X 2 

1 





0 



for t. While this can be solved numerically, a useful 
analytic solution for T g has not been found. However, it 
was discovered that the time at which the system hazard 
rate reaches its maximum is approximately related to the 
time x at which P£(t) reaches its maximum by 



T s « 2x . 



From the numerical results, it appears that T is 

X 2 

related to x by some function of a = •? — , the imbalance 

A 1 

between the two components, i.e. T = g(a) • x. The exact 

s 

form of g(a) is unknown but it appears to be some function 
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whose values are plotted in Figures 2, 3 , and 4 . Note 
that g(a) is never greater than 2.17, nor less than 2.0 
for all 1.001 £ a £ 100, and for values of a > 15, 
g(a) < 2.01. 
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